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ABSTRACT 

The time-averaged 30 ks Chandra AOl HETGS observation of the micro-quasar GRS 1915+105 in the low 
hard state reveals for the first time in this source neutral K absorption edges from iron, silicon, magnesium, and 
sulphur. Ionized resonance (r) absorption from H- and He- like species of Fe (XXV, XXVl) and possibly Ca XX 
are also seen as well as possible emission from neutral Fe Ka and ionized Fe XXV (forbidden, or the resonance 
emission component of a P-Cygni profile). We report the tentative detection of the first astrophysical signature of 
XAFS in the photoelectric edge of Si (and possibly Fe and Mg), attributed to material in grains. The large column 
densities measured from the neutral edges reveal anomalous Si and Fe abundances, and illustrate the importance 
of high resolution X-ray measurements for assessing material which surrounds bright X-ray sources, especially if 
depletion onto grains plays a prominent role. Scenarios for which the anomalous abundances can be attributed to 
surrounding cold material associated with GRS 1915+105 and/or that the enrichment may signify either a highly 
unusual supernova/hypernova, or external supernova activity local to the binary are discussed. We attribute the 
ionized features to a hot disk, disk-wind, or corona environment. Based on H- and He-like Fe (XXV, XXVl), we 
estimate constraints on the ionization parameter (log ^ > 4.15), temperature (T > 2.4 x lO*" K), and hydrogen 
equivalent number density (n > lO'^cm"-') for this region. Variability studies with the simultaneous RXTE data 
show that the light curve count rate tracks the changes in the disk blackbody as well as the power-law flux, with the 
most significant variations observed in the former Similar studies of the Chandra spectra show spectral changes 
which also track the behavior of the light curve, and may point to changes in both the ionizing flux and density of 
the absorber. A 3.69 Hz QPO and weak first harmonic is seen in the RXTE data. 

Subject headings: binary; micro-quasars: general; X-ray: general; individual GRS 1915+105 



1. INTRODUCTION 

The Galactic X-ray transient GRS 1915+105 is an extremely 
energetic X-ray binary that has been studied extensively in 
multiple wavebands which include radio, infrared and X-rays 
(e.g. see review by Mirabel & Rodriguez 1999). Since its 
discovery by Granat/WATCH in 1992 (Castro-Tirado et al. 
1992), it has shown repeated flares separated by intervals in 
a "low-hard" state. In the X-ray band, the overall variability of 
GRS 1915+105 is extremely complex, with a luminosity which 
varies between a few x 10^*^ - few x 10^*^ erg s"' , making it one 
of the most energetic objects known in the Galaxy. Many of 
these variations can be understood in terms of instabilities that 
switch the origin of radiant energy between two components: 
the accretion disk and the X-ray power-law (Muno, Morgan, 
& Remillard 1999). Recent observations in the near-infrared 
H & K band with the VLT reveal the donor to be a K-M III 
giant (Greiner et al,. 2001a) with a 33 day binary period and 
a radial velocity curve interpreted to suggest a compact ob- 
ject of ~ 14M0 (Greiner et al, 2001b). This indicates that 
GRS 1915+105 belongs to the class of dynamically established 
black hole binaries. 

One of the most remarkable aspects of this system is that 
X-ray flares are sometimes followed by superluminal ejec- 
tion events seen in radio images (Mirabel & Rodriguez 1994; 



Fender et al. 1999) — the inferred velocity of the knots in the 
bipolar jets is v ^ 0.95c (i.e. a Lorentz factor of 3). Only a few 
sources in the Galaxy (e.g. GRO 1655^0 and XTE J1748-288) 
display such violent ejections, while other sources exhibit inter- 
vals of more constant radio emission associated with a steady 
jet (Mirabel & Rodriguez 1999). Both types of jet sources qual- 
itatively resemble radio-loud active galactic nuclei (AGN) and 
are often called 'micro-quasars'. 

A key piece to the puzzle for understanding the physical en- 
vironment of GRS 1915+105 is X-ray spectroscopy. ASCA ob- 
servations of GRS 1915+105 have revealed striking structure 
within the (medium resolution) X-ray spectra. Based on obser- 
vations in 1994-1995, Kotani et al. (2000) identified resonant 
absorption lines of Ca XX Ka, Fe (XXV, XXVl) Ka, as weU as 
blends of absorption lines of Ni XXVll Ka + Fe XXV K/? and 
Ni XXVlll Ka + Fe XXVl K/3. However, ASCA is fundamen- 
tally limited in the extent to which it can clearly resolve ab- 
sorption lines or detect emission components. Sufficient spec- 
tral resolution and energy coverage (such as presently possible 
with the Chandra high energy transmission grating) is partic- 
ularly important for obtaining direct information on the line- 
of-sight velocity and velocity dispersion of the material. The 
absolute and relative equivalent widths of the lines can be used 
to determine the column density and ionization parameter of 
the accretion disk environment, and any disk wind. In this 
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paper, we present primarily the prominent spectral features of 
the time averaged X-ray spectrum of GRS 1915+105 obtained 
with the high resolution gratings on-board Chandra. We as- 
sess the chemical abundances based on direct measurements of 
neutral K edges of Fe, S, Si, and Mg (§4). Measurements of 
ionized Fe features are used to assess the conditions of the ac- 
cretion disk corona (§5). Variability is discussed in the context 
of Chandra and RXTE observations (§6). Preliminary results 
were presented in Lee et al. (2001a). 

2. OBSERVATIONS 

GRS 1915+105 was observed with the Chandra High En- 
ergy Transmission Grating (HETGS; Canizares et al., 2001, in 
preparation) on 2000 April 24 beginning at 01:35:49 (MJD: 
51658.06654), with a total integration time ~ 31.4 ks. This 
corresponds to orbital phase ~ 0.8 according to the ephemeris 
quoted by Greiner (2001b), where phase zero (2000 May 02) 
is the time consistent with the blue-to-red crossing in the ra- 
dial velocity. Simultaneous RXTE observations were performed 
with both the Proportional Counter Array (PCA, Jahoda et al., 
1996) and High-Energy X-ray Timing Experiment (HEXTE, 
Rothschild et al., 1998) instruments. 

Observations with the RXTE All-Sky Monitor over the time 
interval of 2000 April 21-28 indicate a relatively steady X-ray 
flux (16% rms variations) and a hard spectrum, with an aver- 
age intensity of 0.4 Crab at 2-12 keV. Radio observations of 
GRS 1915+105 with the Greenbank Interferometer on 2000 
April 24.54 indicate a flux of 20 ± 4 mJy at 2.25 GHz. The 
combined radio and X-ray properties place GRS 1915+105 in 
the "low-hard" state that is associated with the presence of a 
steady jet (Dhawan, Mirabel, & Rodriguez 2000; Muno, Remil- 
lard, & Morgan 2001). 

Both the HETGS and the RXTE PCA instruments detect vari- 
ations from GRS 1915+105 in the form of a smooth, 15% dip 
in the count rates that occurs midway through the observations 
(Fig. 4a, top 2 panels). The spectral variations are subtle, and so 
we first consider the average spectral features and then examine 
the dip in greater detail in §6. 

2.1. The Chandra data & reduction 

The Chandra HETGS is made up of the Medium Energy 
Grating (MEG) with a 2.5-31 A (0.4-5 keV) bandpass, and 
High Energy Grating (HEG) with 1.2-14 A(0.9-10keV) band- 
pass. The peak resolving power (E/AE) for the MEG and HEG 
is ~1000. The resolution of the higher orders improves by a 
factor of n for the nth order, but the spectral bandpass and effi- 
ciency are reduced accordingly. 

We observed GRS 1915+105 during a persistent phase of rel- 
atively low flux. Despite this, the observation suffered from 
photon pile-up. This is the phenomenon by which multiple pho- 
tons fall on a single pixel during the same readout period. As a 
result, the on-board event detection algorithm interprets this as 
a single high energy photon event, and sums the multiple pho- 
ton events to make a single pulse height. Pileup has the effect 
of depleting the photon spectrum, thereby decreasing the pho- 
ton counting rate. Corrections for this effect will be discussed 
in §3. 

Due to severe telemetry and photon pileup problems imposed 
by the large count rate of GRS 1915+105, the observation was 
performed using the graded telemetry mode, and a subarray 
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which reduced CCD frametime by a factor of ~ 2. We addition- 
ally used a 'spatial window' to block out the 0th order image 
(this would have been completely piled-up, even during periods 
of low flux for the energetics of GRS 1915+105). However, be- 
cause the 0th order position is crucial for defining an accurate 
wavelength scale, we calculate its position for our observation 
by finding the intersection of the dispersed MEG-HEG spec- 
tra with the 0th order readout trace. The accuracy of such a 
technique for the determination of the 0th order position is ~ 
0.2—0.3 pixels, which translates respectively to a wavelength 
accuracy of 0.002 A and 0.004 A for flie HEG and MEG 1st 
orders. 

2.2. The RXTE data & reduction 

We have analyzed the RXTE observations of 2000 April 24 
that overlap the Chandra exposure of GRS 1915+105. At high 
time resolution, the PCA telemetry modes yielded 2 ms res- 
olution for two energy channels that span 2-9 keV. We com- 
bined these channels and computed a power density spectrum, 
corrected for instrument deadtime and counting statistics (see 
Morgan, Remillard, & Greiner 1997). We then normalized the 
power spectrum to units of (rms fluctuations/mean)^ Hz"' in the 
source, applying a small correction for the contribution of the 
diffuse X-ray background between 2-9 keV. The results show a 
QPO and a broadband power continuum typical of the low-hard 
state in GRS 1915+105 (e.g., Muno et al., 2001). The strongest 
QPO feature is at 3.69 Hz with an integrated ampUtude of 8.5% 
and a weak first harmonic. 

PCA and HEXTE spectra were derived from standard data 
modes that produce 128-channel PCA spectra (fuU bandwidth) 
every 16 s, and 64-channel HEXTE spectra (full bandwidth) at 
the instrument nodding period of 32 s. There were 3 PCUs in 
operation (units 0, 2, and 3) during the PCA observations, and 
our spectral analyses focus on results from PCUs and 2 be- 
cause of the superior performance of the associated calibration 
files (ftools 5.6), as judged from the analysis of PCA obser- 
vations of the Crab Nebula during late 2000 and early 2001. 

3. THE SPECTRUM 

The best fit continuum model which describes the RXTE 
data is a disk-blackbody plus a power-law component each 
modified by a hydrogen column of 5 x 10^^ cm"^ . The best fit 
color temperature of the disk is 1.4 keV (^ 1.6 x 10^ K) with 
a normahzation of A'bb = 59. (We note that the temperature 
deduced from this black body fit and that in the subsequent 
calculation of §5 based on the ratio of Fe XXV : Fe XXVI are 
not necessarily in conflict given that the DISKBB model is ex- 
pected to overestimate the temperature by as much as a factor 
of 1.7 [Shimura & Takahara 1995], and the value given in §5 
is a lower Umit; however, the temperature measurements here 
and in §5 can also be probing independent components of the 
gas.) The average power-law component exhibits a photon in- 
dex of r ^ 2.4 and a normalization of the power-law at 1 keV 
is 11.4 phcm"^ s"' keV~'. The corresponding apparent (since 
there are losses due to low energy absorption) total 2-25 keV 
flux is 1.89 X 10"^ ergcm^^s"'. The implied (unabsorbed) 
bolometric fluxes are respectively ~ 5 .3 x 10"^ erg cm"^ s"' and 
~ 3 .4 X 1 0~^ erg cm"^ s~^ f or the disk-blackbody and power-law 
(integrated from 1 keV to oo) components, indicating that the 
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power-law component was the dominant emission mechanism 
during these observations. The absolute absorption corrected 
luminosity corresponding to the epoch of our observation is 
lower limit of Ltoi ^ Lx ~ 6.4 x 10^*^ ergs"' at 12.5 kpc. We 
also assess the higher energy data afforded by RXTE HEXTE 
observations, and find that it is best fit with a broken power-law 
with the cutoff fibieak = 30 keV. The photon index below 30 keV 
is r ~ 2.4 (consistent with PC A measurements), and a steeper 
r ^ 2.9 spectrum above 30 keV. Since RXTE is better suited 
for constraining the continuum shape due to its wide band cov- 
erage, we apply this model to the Chandra data. 
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Fig. 1 . — Best-fit Chandra continuum (black) based on best-fit RXTE con- 
tinuum model over-plotted on the (a) non-pileup con'ected Chandra HEG 1 st 
order (blue), MEG 1st order (red), and MEG 3rd (red) order and (b) pileup 
corrected HEG 1st order (blue). 

The overall shape of the continuum discussed above de- 
scribes the Chandra data well (with the exception that the 
best fit power-law and power-law normalization are respec- 
tively r 2.1 and 5.3 phcm"^ s"' keV~' for fits to the Chan- 
dra data), especially given the problems with pileup, which is 
responsible for ^ 40% of the 60% normalization difference 
between the RXTE and Chandra data. (Pileup will affect differ- 
ent parts of the spectrum by different amounts.) To illustrate. 
Fig. la shows the absorbed disk-blackbody plus power-law 
model over-plotted on the HEG 1 st and MEG 1 st & 3rd order 
data. Spectral deviations from the model in the MEG/HEG 1st 
order spectra at « 2-5 A is due to pileup as a result of the high 
count rate of the source in this passband; instrumental pileup 
severely affects the 1st orders between « 6.5-9 A since this is 
where the effective area of the HETGS peaks. The MEG 3rd 
order spectrum at < 2 A is affected by both first order photon 
pileup at ~ 5 A (i.e. MEG 3rd order pileup = A/order = 5A/3), 
and 0th order scattering. We extract a separate set of HEG and 

^http://space.mit.edu/CXC/ISIS/ 
'http://www-cxro.lbl.gov/optical_constants/asf.html 



MEG 1 st order spectra to be used with the HETG pileup model 
as implemented in ISIS^ (see Davis 2001, Houck & DeNicola 
2000). These data are extracted with a coarser (respectively 
for the HEG and MEG 1st orders -0.017A and 0.033A) bin- 
ning to reflect the size of an ACIS event detection cell. The 
pileup corrected data shows a reduction from 60% to a few % 
in the overall normalization offset between RXTE and Chandra 
(Fig. lb). We note that the broadband RXTE PCA flux is typ- 
ically ~ 15-25% higher when compared against other X-ray 
missions (e.g. ASCA , HEXTE). 







Mg K Edge 


MEG 1st order 








Wavelength (A) 

Fig. 2. — Photoelectric K-shell edges of S, Si and Mg. Overplotted are the 
respective best-fit models described in §2. Note the prominence of the Si edge, 
and possible XAFS structure. (The Fe K edge is shown in Fig. 3.) 



4. PHOTOELECTRIC EDGES 

We find evidence for prominent absorption edges due to K- 
shell Fe, S, Si, and Mg (Table 1, Figs. 2 & 3). Because the spec- 
trum of GRS 1915-1-105 is complicated by ionized resonance 
features (e.g. Fig. 3), and possibly resonance structure due to 
X-ray absorption fine structure (XAFS), we derive t by fitting a 
linear function to pieces of the data bracketing ~ 0. 1 A to either 
side of a given edge, and extrapolate the fit to the expected edge 
energy (see Fig. 2). (XAFS in the spectrum of GRS 1915H-105 
are discussed later in this section.) Additionally, since pileup 
will also affect the depth of an edge by reducing the jump in the 
discontinuity at the edge, we use the pileup corrected 1st order 
HEG to assess the optical depth of Fe, S, & Si, and MEG to 
assess Mg. The column densities (shown in Table 1) are calcu- 
lated using the value for the K-shell photoionization cross sec- 
tion of the relevant species ( ^ Henke et al., 1993), and evaluated 
for solar and ISM abundances (Wilms, Allen, McCray 2000). 
Given the possibility for emission and absorption features and 
edge structure which we have not accounted for, a conservative 
estimate for the systematic uncertainties in r is calculated us- 
ing the error in the counts and best fit continuum level on either 
side of the edge to obtain a fractional uncertainty. 
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Table 1 



Prominent Neutral K Edges of the Time- averaged data 







Photoelectric Edges 




Species 


Ao 

(A) 


T 


"Nh (solar) 
(10"2 cm-2) 


''Nh (ISM) 
(10^2 


MgK 
SiK 


9.52 
6.74 


0.29 ± 0.03 
0.47 ±0.01 


3.1+^-^ 
3 4^:? 


-OA 

16 l"*^ ^ 
10.1 „ 3 

4 fi^-^ 
^•°-o.9 

10.9!!:^ 


SK 
FeK 


5.02 
1.74 


07+0.002 
"•"'-0.01 
0.11 ±0.02 


3:2^8:? 

Q ■xA'^ 
^■^-13 



"'^ The total hydrogen column assuming (a) solar and (b) ISM 
abundances. The expected atomic (Wh) and molecular (Nm) hydro- 
gen along the line of sight is respectively > 1.76 x lO^^cm"^ and 
~ 7.43 X 10^' cm-^ . The MEG summed 1st order data is used to assess 
the Mg edge; all other edges were fit using the HEG summed 1st order 
data. 



We compare our results (below) with constraints on the 
Galactic A'hg ^ 1-8 x 10^^ cm"^ (determined from 21 cm 
emission from atomic hydrogen along the line of sight, Dickey 
& Lockman, 1990). Despite the corrections for optical depth 
effects, this value is regarded as a lower limit because of its 
reduced sensitivity to both cold and ionized H, as well as H2. 
Dame et al. (2001) report additional molecular hydrogen (based 
on CO maps) withA^H2 ~7.43 x 10^' cm"^ to the nearest to- 
wards GRS 1915H-105. Therefore, the line of sight column can 
be > 3 X 10^^ cm"^ with contributions from both atomic and 
molecular hydrogen. (Chapuis and Corbel 2001, in prepara- 
tion, find the total [atomic plus molecular] line of sight column 
density toward GRS 1915H-105 from millimeter and HI obser- 
vations to be 3.3 ± 0.2 x 10^^ cm"^ ; private communication S. 
Corbel.) 

Using the HETG absorption edges for S and Mg, we derive 
A'h values ~ 3 x 10^^ cm~^, assuming solar abundances, and 
^ 5 X 10^^ cm"^, assuming ISM abundances. These values are 
consistent with the A^h value derived from low energy contin- 
uum absorption. However, the A^h derived from Si & Fe (for 
solar abundances) are greater by > 5 x 10^^ cm~^ (and much 
greater for ISM abundances) when compared with what is ex- 
pected from the line of sight column density (Table 1). We note 
that earlier reports by Lee et al. (2001a) for A'h were based on 
data that had not been corrected for pileup. 

We point out the exciting possibility that astrophysical XAFS 
are present in the spectrum of GRS 1915-1-105, most noticeable 
in Si, but may also be present at the Mg and Fe edges. XAFS is 
a general term for both EXAFS (extended X-ray absorption fine 
structure) and XANES (X-ray absorption near edge structure), 
with the distinction based on single (EXAFS) versus multiple 
(XANES) scatterings of the photoelectrons from the immedi- 
ate neighboring atoms (e.g. Sayers, Stern, & Lytle 1971, Woo, 
Forrey, & Cho 1997, and references therein). Such features 
appear as wave-like modulations, with the former manifest in 
structure at > 50 eV from the absorption edge, and the latter 
in the region between the edge and EXAFS region. A quali- 
tative comparison of Si in the detector (Si02 and polysilicon) 
shows that the modulations of the Si EXAFS in our observation 
are much stronger than that expected from the detectors, and 
therefore may be features which are associated with the X-ray 
source. We note further that the detector Si region has been 
carefully calibrated against ground based measurements at syn- 
chrotron beamlines (e.g. Prigozhin et al., 1998; also Mori et 

"'http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 

"CXC Proposers' Observatory Guide, December 2000, Version 3, p.l90 



al., 2001 for further studies of absorption edges in front illu- 
minated CCDs). Therefore, we believe that a definitive claim 
of XAFS in GRS 1915-1-105 is hampered more by the lack of 
adequate statistics than calibration errors in the quantum effi- 
ciency of the Chandra CCDs. (The calibration near the Si edge 
is good to ~ 3% over 0.02 A bins - private communication, H. 
Marshall.) The possible EXAFS at the Si edge are noted in the 
instrument-corrected spectrum shown in Fig. 2. An investiga- 
tion of several slightly less absorbed X-ray binaries which have 
good statistics and for which the instrumental response was ap- 
plied (Schulz et al., 2000) do not show XAFS such as that seen 
here, and gives us further confidence that the XAFS seen in the 
Chandra spectrum of GRS 1915H-105 is associated with an as- 
trophysical phenomenon. At present, we merely point out this 
exciting possibiUty - an in-depth treatment of XAFS is beyond 
the scope of this paper, and will be addressed in a subsequent 
paper which will investigate the presence of XAFS in highly 
absorbed bright X-ray binaries. 

We further verify that the Si edge and possible structure is 
intrinsic to the source by comparing the measured edge depth 
against the contribution from the Chandra CCDs. The disconti- 
nuity in the effective area at the detector Si edge is ^ 20-30%, 
and verified to be accurately " modeled to < 10% but as stated 
previously can be good to ~ 3% over 0.02 A bins; the latter is 
particularly important for the assurance that the residuals which 
we attribute to XAFS at the Si edge are not likely to be an in- 
strumental feature. This is compared with the > factor 2 dis- 
continuity measured from the HEG 1st order count spectra. We 
find further confidence that the Si edge is not a detector fea- 
ture from the MEG 3rd order spectrum, where the edge remains 
prominent in the data, but for which there is a minute contribu- 
tion from the detector. (The spectra for the relevant orders are 
the combined plus and minus sides; for the MEG 3rd order, the 
Si edge seen in our spectrum for one of the sides fall on the 
backside (SI) chip where there is no/little contribution from the 
detector Si edge. - See e.g. the effective area curves shown in 
Fig.2ofLeeetal.,2001a) 

5. THE IRON COMPLEX - A HOT DISK ATMOSPHERE ? 

In addition to neutral edges, and possible neutral (Fe Ka) 
emission, the spectrum of GRS 1915H-105 also exhibits a strong 
Fe XXVI resonant absorption feature (Fig. 3). A weak Fe xxv 
resonant absorption may also be present. We detect a ^2.1 a 
emission feature at ~ 1.87 A which can be either the Fe xxv 
Is^^Sq - l5'25'-'5'i (forbidden) emission, or the emission com- 
ponent of a P-Cygni resonance line of Fe XXV (similar to what 
is seen in Cir X-1 - Schulz & Brandt, 2001, in preparation). Be- 
cause the Fe region is < 5% piled-up, we use the original HEG 
combined plus and minus 1st order data binned to 0.0075 A to 
assess the features detailed in Table 2. The ~ 0.002 A wave- 
length errors are consistent with the expected accuracy for the 
HEG data (e.g. §2.1), and present measurements show consis- 
tency with a zero velocity shift. (The ASCA spectrum of the 
microquasar GRO 1655-40 (Ueda et al., 1998) also shows no 
evidence for blueshifted lines.) The resolution of the Chandra 
HEG is ~ 1600 kms^' FWHM (-10 mA) at the iron energies. 
The velocity widths of Table 2 indicate that the lines are ei- 
ther not resolved or only marginally so. The larger breadth to 
the Une width of the Fe xxvi absorption when compared with 
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Prominent Absorption & Emission Features in the Time- averaged data 



Species 


" A 




'AA 


''a 


Flux 


' ' re 






(A) 


( km s"' ) 


(kms"' ) 






A C /V /V V 1 -i— 'V CI. 


1 77Q8 


1 j-j^j +0 0074 

i. . / / _I_ Vj . VJVjZ-T^ 


_770 + 400 


'C7Q +400 


-1.72 ±0.41 


4.7 


Fe XXV Ly a*' 


L8505 


L8547± 0.0062 


+680 ± 1000 


= Fe XXVI Ly a 


-0.61 ±0.38 


1.7 


Fe XXV Ly a 


L8505 


L8540± 0.0070 


+570 ±1140 


499 ± 728 


-0.58 ±0.37 


1.6 


Fe XXV 


L8682 


L8676± 0.0027 


-100 ±430 


= Fe XXVI Ly a 


0.98 ±0.36 


2.7 


Fe XXV t 


L8682 


L8685± 0.0028 


+50 ±450 


164 ±346 


1.00 ±0.36 


2.7 


FeKa 


L9361 


L9350± 0.0094 


-170 ±1457 


380 ±670 


-0.71 ±0.35 


2.0 


Un-IDed 




L7315 ±0.0016 




177 ±165 


< 1.96±0.51 


<6 


Ca XX Ly a*^ 


3.0203 


3.01 17 ±0.0029 


-850 ±290 


205 ± 135 


>0.34±0.17 


> 1.2 



Table 2 

'Absorption feature; otherwise emission. ^Line is identified as either the Fe XXV (/) emission or the emission portion of the Fe XXV (r) 
P-Cygni. 'The line width is tied to that of Fe XXVI absorption. ^Lower limit values quoted since region is heavily piled-up. " Rest wavelengths of 

identified hnes from Verner et al., 1995 ^ Best measured centroid wavelengths (A). Amount of shift of measured wavelength from rest (- for 
blueshift, + for redshift) 'True' velocity width (km/s). The instrumental line spread function has already been subtracted out. Flux in the line. 

' Equivalent width. 



the A ~ 1.87A Fe XXV emission (if the forbidden line) may 
indicate that the region responsible for the Fe XXV emission 
may only be part of the region responsible for the Fe XXVI ab- 
sorption. Alternatively, the A ^ 1.87A line is the (r) emission 
component of the Fe XXV P-Cygni complex which is equally 
divided between absorption and emission. 

Because the Fe XXVI and Fe XXV resonant absorption fea- 
tures are detected, we can estimate the physical parameters of 
the plasma from which they originate based on the ratio of the 
Fe XXVI : Fe XXV column densities (A^Fe26 : A^Feis)- If the lines 
are resolved (e.g. Fe XXVl), we can assume that they are unsat- 
urated and on the linear part of the curve of growth, and there- 
fore estimate the Fe XXVI and Fe XXV column densities to be 
A^Fe26 4 X 10'^ cm-2 and A?Fe25 7 X 10'^ cm based on the 
relation 



^ = ^f'i = 8-85 X lO-'^A^,. Xfij (1) 



Q I I I I I I I I I I I I I I I I I I I I I I I I ' I I 

Fe XXV -> 




Fe K Fdge 
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Fig. 3. — Time-averaged AOl 30ks HEG 1st order data of the Fe region 
binned to 0.0075 Ashow both ionized (most likely from the accretion disk at- 
mosphere) and neutral (from the source environment, or along the line of sight) 
Fe features. Over-plotted is the best fit continuum (dashed blue) and identified 
(red) lines. The unidentified line may be a shift in the edge due to XAFS. 



where Nj is the column density for the relevant species, fij (re- 
spectively for Fe (XXV, XXVI) is 0.798, 0.416, from Verner 
1995) is the oscillator strength, W\ is the equivalent width of 
the line, and A is the wavelength in cm units (Spitzer 1978). 
(If instead, the features are due to a superposition of satu- 
rated lines over a range of velocities, the column densities 
quoted would be higher) The Icr confidence upper limit for 
the Fe XXV flux in Table 2 gives an equivalent width ^2.7 
mA , and therefore A^Fe25 ^ 11 x 10'^ cm"' (compared with the 
bestfit A^Fe25 7 X 10"' cm"^ quoted above). We use this upper 
limit value of A^Fe25 for our subsequent calculations to estimate 
the physical conditions (e.g. temperature, ionization parameter, 
etc.) of the hot plasma. Results based on calculations using the 
best-fit Fe XXV flux and equivalent width will be noted in '[]'. 



Based on the calculations of Kallman & Bautista (2001) for 
a power law ionizing spectrum of energy index a = -l (F = 2), 
and photoionized plasma, the ratio of A^Fe26 to A^Fe25 implies 
an ionization fraction of < 0.45 0.39] for Fe XXVI from 
which we estimate log ^ > 4.15[ergcms~'] [log ^ ^ 4.3], and 
T > 2.4 X 10^ K [- 2.7 x 10^ K]. (This is a much higher ion- 
ization state than that reported by Kotani et al., 2000 for their 
high temperature limit.) Given the likelihood that these ion- 
ized features originate in close proximity to the bright X-ray 
source, we expect photoionization (rather than collisional ion- 
ization) to be the dominant physical process, due to the high 
radiation field. (Temperatures of > 10*^ K would be neces- 
sary for producing hydrogen-like iron in a coUisionally ion- 
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ized plasma.) We next correct for the ionization fraction to 
obtain an equivalent hydrogen column A^h in this region re- 
spectively for solar and ISM abundances to be 2.8 x 10^^ cm~^ 
and 3.3 X 10^^ cm"-^ . Based on the relation ^ = L/n/?^ erg cms"' 
where A^h = nAR = the maximum column density in this re- 
gion, we estimate R^/AR ~ 2 x 10^^ cm. (R is defined to be 
the distance from the ionizing flux to the absorber and AR the 
thickness of the absorber.) If we further postulate that the vol- 
ume filling factor AR/R must be small (e.g. « 0.1) in order 
that ^ not change over the region, we estimate an upper limit on 
/? < 2 X 10" cm, implying that A/? « 2 x 10"' cm. It follows 
that the hydrogen equivalent number density n > 2 x lO'^cm"^ 
(assuming solar abundances). 

It is plausible that the ionized absorber responsible for the 
H- and He-like Fe lines at a distance R < lO" cm is associated 
with material flowing out from the X-ray source (e.g. wind, see 
also §6), which we can compare with the accretion rate. Us- 
ing the relation for the ionization parameter, we can define the 
(spherical) mass outflow rate M/iow at a velocity v ~ 10^ cm s~' 
(100 kms-i)tobe 

Mfiow = ^T^r-nmpV{—) = 4nmpV(—)( — ) (2) 
47r 4 ^'^ 

~9.5 x 10'^(-^)gms-' 

where r is the characteristic radius, p = nnip is the density of 
the absorbing material, with n equal to the number density of 
electrons in this material and Wp the proton mass, and fl = solid 
angle subtended by the outflow. This can be contrasted with 

Maccretion = ^ ~ 7. 1 X lO'^ gms"' (3) 

where the efficiency 77 0.1 and the lower limit on the bolo- 
metric luminosity Lhoi ^ is used, (from §3). Such a compar- 
ison for V ^ 100km s"' shows that as the covering fraction (i.e. 
r2/47r) approaches unity, the mass flux in the flow (or wind) can 
be comparable to that in the disk. 

The high ionization state of the hot medium (even based on 
the lower limit log ^ ^ 4.15 ) which supports the H- and He- 
Uke Fe lines is too ionized to support most other species with 
the exception of H-hke Ar xviii (ionization fraction ~ 9%) and 
H-Uke Ca xx (~ 16%), and Fe xxv (~ 12%), compared to the 
^A5% ionization fraction for Fe XXVI. (These percentages will 
decrease with increasing ionization.) The lower limit (since this 
region is severely affected by pileup) to the column density of 
Ca XX (Jij = 0.416) is A^caio ^ 3.6 x 10'* cm'^ 

6. TRACKING VARIABILITY IN THE DISK ? 

The count rates from the Chandra HETGS and the RXTE 
PCA instruments for GRS 1915-1-105 are shown in the top two 

panels of Fig. 4a. To investigate the physical processes which 
may be related to the observed temporal variability, we use the 
PCA spectral parameters to compute the inferred bolometric 
flux from the accretion disk (Fig. 4a, panel 3) and from the 
power law component (1-20 keV; panel 4). We find that the 
most significant variations occur in the disk blackbody compo- 
nent, although there is also a smaller dip in the power-law flux 
as well. 

The time-sliced Chandra data shown in Fig. 4b suggests that 
evolution of the ionized line features may also be correlated 
with changes to the disk flux. In particular, the Chandra spec- 
tra shows that the Fe xxv line is more prominent during the dip 
between the ~ 10-20 ks of our observation, in contrast to the 



prominence of the Fe XXVI absorption at the beginning and end 
of the observation when the source is brighter. The time interval 
corresponding to the strongest Fe XXV absorption is coincident 
with the dip in the Fig. 4a (panels 1 & 2) light curve count 
rate, as well as the notable diminution of the disk flux (panel 
3), and to a lesser extent, the power law flux (panel 4) seen in 
the RXTE data. This would support the idea that the ionized 
absorption features originate from the disk. However, we note 
that an ~ 20% change in the flux at the ionization discussed 
in §5 implies an ~ 30% change in the ionization fraction of 
Fe XXV, if we assume that ^ scales with the ionizing flux, while 
n and R remain constant. This is not sufficient to explain the 
lack of Fe XXV absorption, and dominance of Fe XXVI absorp- 
tion during the periods when the count rate is high. An ~ 50% 
change in ^ which will reduce the ionizing fraction of Fe XXV 
by a factor of two might be more plausible for explaining what 
is observed, indicating that the ionizing flux is only part of the 
solution, and changing density the other factor. If we consider 
the calculations in §5 for a wind with velocity v ~ 10^ cm s~' 
and distance to the absorber/? ~ lO" cm, it is plausible that the 
spectral variabiUty can be attributed to a flow which can change 
on timescales of 10 ks such as that seen here. This is not in- 
conceivable since 10 ks is a very long time, corresponding to 
many dynamical timescales of the relevant parts of the disk, in 
the life of a source Uke GRS 1915-H105. Furthermore, since this 
source accretes near its Eddington Umit, one would expect the 
mass flow rate and wind density to be strong functions of the 
luminosity, such that the actual structure of the wind (i.e. den- 
sity and size) can change as a function of the luminosity. Simi- 
larly observed variations in Cir X-1 have been interpreted to be 
due to changes in the ionization fraction of the wind (Schulz & 
Brandt, 2001, in preparation). We note however that the statis- 
tics are not adequate for us to present a complete quantitative 
picture for GRS 1915H-105, especially in light of the numer- 
ous problems we have with pileup that need to be considered. 
Cycle 2 Chandra data of GRS 1915H-105 taken in a mode that 
optimally reduces pileup will be able to study the variability 
aspects in greater detail. 

7. DISCUSSION 

The Chandra HETGS spectrum of GRS 1915H-105 shows 
neutral and ionized features superposed on a continuum best 
described by a disk blackbody and a dominating power-law 
modified by a column of ~ 5 x lO^^cm"^ absorption. The X- 
ray spectral and temporal characteristics as well as the 3.7 Hz 
QPO seen in our data, and the detection of faint radio emis- 
sion are all consistent with the low-hard state. This indicates 
that the X-ray power-law is the dominant emission component 
and there is a weak and persistent radio jet. No obvious red- / 
blue-shifted features indicative of a jet can be seen in the time- 
averaged Chandra spectrum. 

The H column density derived from the Si and Fe edges 
(assuming solar abundances) is ~ 8-9 x lO'^'^cm"^, which is 
substantially higher than the values derived from S and Mg, 
and continuum absorption determined from broad-band instru- 
ments. The discrepancy for excess Si and Fe absorption is even 
higher if we assume ISM abundances. 

There is a possibility that the high X-ray column density and 
abundance excesses may be related to material that is associated 
with the immediate environment of the source. Both Mirabel et 
al. (1996) and Marti et al. (2000) have suggested the possibil- 
ity for heated dust that is associated with the surroundings of 
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Fig. 4. — (a) Simultaneous Chandra (top panel 1) and RXTE (panel 2) light curves during the epoch beginning 2000 April 24 at 01:35:49 (MJD: 51658.06654). 
Changes in the blackbody flux (panel 3) and power-law flux (panel 4) are coincident with changes in the light curve count rates. (Panels 3 & 4 are in units of 
10"** ergcm"^ s"' .) (b) The 30 ks Chandra spectrum divided into three 10 ks intervals depicts the evolution of the spectral features. Note in particular the promi- 
nence of the Fe XXVI absoiption at the beginning and end of the observation when the flux was highest, and of Fe XXV during the dip. Over-plotted is the rough 
shape of the continuum (red-dashed) from the time-averaged spectra shown in Fig. 3. Error bars are representative of a statistically independent bin for which the 
data have been boxcar smoothed by 3, coiresponding to ~ 0.0075 A. 



GRS 1915-1-105 from IR observations. If true, these anomalous 
Si and Fe abundances may provide hints to the evolutionary 
events which gave rise to this source, and the black hole in par- 
ticular. However, we acknowledge that more data from other 
X-ray binaries with similar galactic column densities are nec- 
essary to determine the absorption abundance ratios over the 
Galaxy and so confirm that the results for GRS 1915-1-105 are 
exceptional. Thus far, observations of SS 433 (Marshall et al. 
2001) and Cyg X-1 (Schulz et al. 2001) do not show such Fe 
and Si overabundances. (See also e.g. respectively Paerels et al. 
2001 and Schulz et al. 2001 for the X-ray binaries X0614H-091 
and Cyg X-1, and Lee et al. 2001b for the extragalactic source 
MCG-6-30-15.) Clearly, further sensitive measurements such 
as these of photoelectric edges with Chandra will provide valu- 
able information for assessing the material which surrounds 
bright X-ray sources especially if depletion plays a prominent 
role. Such measurements will be particularly important for con- 
firming the tentative detection of XAFS (particularly in Si) as 
reported here. If true, such a positive detection of these absorp- 
tion fine structure representative of material in solids (as e.g. 
interstellar grains) will have important consequences for con- 
ducting solid state astrophysics, by which grain properties can 
be extracted via the solid's inner compound structure. 

In the interim, the abundance excesses for X-ray absorption 
edges in GRS 1915-1-105 can be contrasted with the discovery of 
abundance anomalies in the optical spectra of companion stars 
in the microquasars GRO J1655-40 and V4641 Sgr (Israelian et 
al. 1999; Brown et al. 2000; Orosz et al. 2001). In the latter 
cases, the over-abundances of a-process elements, but not the 
Fe group, has led to the hypothesis that the companion star cap- 
tured supernova or hypernova ejecta is related to the formation 



of the black hole. The Fe in the progenitor's core is likely to 
be efficiently captured by the black hole. In the case of GRS 
1915-1-105, the Fe and Si overabundance pertains to material 
along the line of sight, and the enrichment may signify either a 
highly unusual supernova or perhaps external supernova activ- 
ity local to the binary. 

One way to reconcile the overabundance of Fe in 
GRS 1915-1-105 is if an asymmetric hypernovae/supernova ex- 
plosion produced heavier elements like Fe and Si along the 
poles. For the UMq black hole in GRS 1915H-105, the binary 
spin axis would then have to be misaligned with the 70° in- 
clination of the jet (Mirabel Sl Rodriguez 1994) with respect 
to the line of sight. Such a misalignment could be induced by 
disk warping, precession of the spin axes, or residual angular 
momentum from the initial supernova explosion. 

One rather speculative possibility is that the excess Fe and Si 
might indicate a connection between some kind of SN and BH 
microquasars. Given that BH microquasars are likely formed 
from SN, Gamma-ray bursts (GRB) may also represent some 
kind of microquasar (e.g. Blackman et al. 1996; Kouveliotou 
et al. 1999) formed similarly (Paczynski 1998; Brown et al. 
2000). This makes contact with the results of Brown et al. 
(2000) for GRO 1655H-40, although as mentioned above, there 
the overabundant elements are alpha-process elements S and O 
which are more easily understandable as a SN signature. Nev- 
ertheless, one might ultimately imagine a classification of mi- 
croquasars analogous to that of AGN. 

In addition to the first detection of photoelectric edges in 
GRS 1915-1-105, our Chandra observation confirms the pres- 
ence of absorption features from highly ionized gas previously 
reported by Kotani et al. (2000) with ASCA data taken between 
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in 1994 to 1999. Despite their use of a very complicated con- 
tinuum model, it is evident from Fig. 1 of Kotani et al. (2000) 
that the observations showing strong absorption features (years 
1994 and 1995) display steep X-ray spectra near 10 keV, in con- 
trast with the 1996 observation, where the spectrum is harder, 
and no absorption features were found with the ASCA data. We 
confirm with the RXTE PCA data archive that GRS1915+105 
was in a bright mode (1.2 Crab) of the hard-steady state on 1996 
October 23 (see Muno et al. 2001). One might therefore suspect 
from the reported ASCA observations that the absorption fea- 
tures from highly ionized atoms are confined to the soft X-ray 
states of GRS 1915+105; however, the Chandra results shown 
here clearly show such absorption features in the low-hard state 
(0.4 Crab), as well. This is consistent with the view that the ab- 
sorbing gas is associated with the accretion disk, which appears 
evident as the soft component in the X-ray spectrum. 

Given that GRS 1915+105 is probably accreting near or 
above its Eddington limit, the presence of a strong radiatively- 
driven wind is almost expected. However, we do not find 
obvious P-Cygni profiles attributed to a disk wind — i.e., 
redshifted/rest-frame emission lines (from material out of the 
line of sight) accompanied by blueshifted absorption lines 
(from the foreground, line-of-sight parts of the wind) in the 
time-averaged data, akin to what has been reported for Cir X- 
1 (Brandt & Schulz 2000). The lack of a wind from the star 
may corroborate the finding that the donor in GRS 1915+105 
is a LMXB (Greiner et al., 2001a) where accretion via Roche 
lobe overflow is the most plausible mechanism for mass trans- 
fer. However, the wind (and/or jet signatures) if present (e.g. 
§5) in GRS 1915+105 is Ukely to originate from the disk, and 
may show up in a more subtle form (e.g. through variability 
studies, §6). The ionized Fe absorption lines of the time aver- 
aged spectrum are likely to have origins from the disk or disk 
atmosphere. 

This hot absorbing medium is a particularly important topic, 
since it points to a disk atmosphere, and possibly even disk 
wind, or some other type of hot gas structure that is absent 
from most models for black hole accretion. It will be inter- 
esting to diagnose the effects of scattering and absorption in 



this medium with Chandra spectra of GRS 1915+105 in dif- 
ferent flux states, in order to build a picture of which spec- 
tral changes are inherent to the source and which are tied to 
changes along the line of sight to the source. The latter topic 
is a highly controversial one for GRS 1915+105. The violent 
changes seen in the apparent disk radius and temperature have 
been interpreted on some occasions as physical changes that 
represent the disk instability (Belloni et al 1997) and the ef- 
fects of the jet production mechanism (Eikenberry et al. 1998). 
However, analogous disk changes, particularly when the vari- 
ations are more modest, have been alternatively interpreted as 
effects due to radiative transfer in an accretion disk (Merloni 
Fabian & Ross 2000). The Chandra spectra (§6) suggest that 
the apparent changes in the strength of the Fe XXV and Fe XXVI 
absorption can be linked to changes in the incident ionizing flux 
and density, and may point to a flow which can change over 
~ 10 ks time scales. Clearly, longer continuous observations 
with Chandra and RXTE are needed for this source in order to 
monitor the observed variability, and differentiate between the 
different mechanisms which may be responsible. 
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